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Abstract 

O |. Wc present an empirical correction of sky coordinates of X-ray photons obtained with the X-ray Imaging 

Spectrometer (XIS) aboard the Suzaku satellite to improve the source position accuracy and restore the 
point-spread function (PSF). The XIS images are known to have an uncertainty in position of up to 
1', and to show considerable degradations of the PSF. These problems are caused by a drifting of the 
satellite attitude due to thermal distortion of the side panel #7, where the attitude control system is 
mounted. Wc found that the position error averaged over a pointing observation can be largely reduced 
by using the relation between the deviation of the source position in the DETX direction and the ecliptic 
latitude of the pointing target. We parameterized the wobbling of the source position synchronized with 
QO ' the 96-minute satellite orbital period with temperatures of onboard radiators and elapsed time since the 

night-day transition of the spacecraft. We developed software, aeattcor, to correct the image drift using 
these parameters, and applied it to 27 point-source images taken in the Suzaku initial science-operation 
phase. We show that the radius of the 90% error circle of the source position was reduced to 19" and the 
' PSF was sharpened. These improvements have enhanced the scientific capability of the Suzaku XIS. 
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LJ' 1. Introduction terparts in other wavelengths [e.g., a new supersoft source 

discovered by Suzaku; see Takei et al. (2007)] and also 

The X-ray Observatory Suzaku (Mitsuda et al. 2007), to construct source catalogs [e.g., Ueda et al. (2001) and 

JH \ developed jointly by Japan and the US, has two ac- Ueda et al. (2005) in the case of ASCA]. 

tive scientific payload, XIS (X-ray Imaging Spectrometer: The celestial coordinates of each photon detected with 

Koyama et al. 2007) and HXD (Hard X-ray Detector: the Suzaku XRT-XIS are calculated from the position 

Takahashi et al. 2007) . XIS consists of four X-ray imaging coordinates of the CCD pixels defined on the XRT fo- 

CCD cameras covering a field-of-view of 17.'8 x 17f8; three cal plane and the instantaneous pointing direction of the 

are front-illuminated (0.4-12 kcV), and one is a back- XRT. The arcsecond localizations of optical guide stars 

illuminated CCD (0.2-12 kcV). The XISs are located on by a pair of star trackers (STTs) primarily determine the 

the focal plane of the X-ray Telescopes, XRTs (Scrlcmitsos absolute pointing direction of the XRT. While guide stars 

et al. 2007). Each X-ray photon recorded by the XIS can arc not available, a set of gyroscopes takes a role to deter- 

be localized with the CCD native pixels, each of which mine the pointing direction. However, over the course of 

has dimensions of 1.04" x 1.04". While the PSF of the in-orbit calibrations of the Suzaku XRT, it was revealed 

XRTs has a broad distribution with a half-power diam- that the calculated positions of X-ray sources have an er- 

eter (HPD) of ~ 2', it has a sharp core with a width of ror of 1' (Serlemitsos et al. 2007), which is significantly 

~ 10". Thus, the sky positions of point-like sources with larger than that expected in the pre-fiight specifications, 

good photon statistics 10 4 photons) are in principle de- ~ 10". The deviation of the calculated source position 

termined with a precision of better than ~ 10", if we know from the true position is time-dependent and synchronized 

the instantaneous pointing direction of the XRT perfectly, with the orbital motion of the spacecraft. Its behavior 

The positional information is crucial to search for coun- differs from one observation to another. The error of the 
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source position can be as large as 1', and therefore the 
astrometric accuracy of the XIS images has been taken to 
be ~ 1' so far. 

The cause of this problem has been identified with 
a thermal distortion of the satellite chassis, particularly 
the cross frames on the side panel #7 (see figure 1), as 
briefly stated in Serlemitsos et al. (2007). This panel 
is situated on the anti-Sun side of the spacecraft, where 
the dewar of the X-ray Spectrometer (XRS: Kelley et al. 
2007) is placed. To expose the XRS dewar to space, the 
lower half of the panel #7 is skeletonized by cross frames 
made of aluminum [see figure 1 and also Mitsuda et al. 
(2007)]. The cross frames suffer from thermal distortion 
due mainly to Earth albedo illumination. Specifically, side 
panel #7 is forced to be tilted in the DETY direction as 
the cross frames expand, since the upper end of the panel 
is fixed at the side panels #6 and #8. Since the STT and 
gyroscopes are co-located on side panel #7, it causes a 
time-dependent misalignment between the XRT-XIS sys- 
tem and the attitude system, and consequently drifts the 
XRT pointing in the DETY direction. Also, an imbalance 
of thermal expansion between the cross frames induces ro- 
tation around the Y-axis (perpendicular to the side panel 
#7), resulting in drifts in the DETX direction. 

We have investigated the characteristics of drift of the 
X-ray images to restore the pointing accuracy as close to 
the pre-flight specification as possible. Suzaku docs not 
have onboard sensors for the temperatures of the cross 
frames, which arc considered to be the most relevant pa- 
rameters to the thermal distortion. We therefore have 
searched for relevant key parameters within the limited 
housekeeping data available on ground. 

In the present work, we parameterize the drifts of the 
XRT pointing direction against the attitude control sys- 
tem with the spacecraft temperatures, the satellite orbital 
phase, and the average pointing direction relative to the 
ecliptic plane, and found a method to correct the position 
accuracy and the PSF in the observed XIS image using 
these parameters. We applied the method to point-source 
images taken in the initial scientific operation phase and 
derived the position accuracy currently achievable. The 
outline of this paper is as follows. We describe the data 
sets and basic reduction in section 2. The parameteriza- 
tion of drifts of the XRT pointing is presented in section 3, 
which is followed by a brief description of new software, 
aeattcor, developed for the position-error correction in 
section 4. We present the improvement on the PSF and 
the position accuracy in section 5. A summary is in sec- 
tion 6. 

2. Observations 

We analyzed XIS data acquired during the Suzaku ini- 
tial science-operation phase for the present work. We 
selected data of point sources whose celestial coordi- 
nates were cataloged with a precision of better than 
a few arcseconds, their X-ray fluxes were higher than 
0.5 counts s -1 chip -1 , but less than 100 counts s _1 chip -1 
because a too-high count-rate caused a significant pileup 
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Fig. 1. Schematic view of the Suzaku satellite. Some key 
components in this work are indicated by arrows. Side panels 
#6 and #8 are indicated by grey shaded areas. The space- 
craft (S/C) XYZ directions are shown. The DETX direction 
coincides with the S/C +X direction, while DETY is defined 
by the S/C — Y direction. 

at the peak of the PSF and degraded the precision of the 
position determination. 

Table 1 lists the selected X-ray sources, consisting 
mainly of active galactic nuclei and X-ray binaries, used 
in this work. (For multi-epoch monitoring observations, 
we used only the data of the two latest epochs in order not 
to enhance statistical weights on these sources.) We re- 
trieved the Suzaku XIS data with version 1.2 pipeline pro- 
cessing (V1.2), and utilized cleaned event data. We ana- 
lyzed the XIS images using sky coordinates, X/Y columns 
(+X is directed to west and +Y is north) in the event 
files, which were calculated from the detector coordinates 
(DETX/Y) of the CCD pixels and the satellite attitude 
monitored by the STTs and gyroscopes (Ishisaki et al. 
2007). It is worth emphasizing that prior to the pipeline 
processing VI. 0, the attitude determination with STTs 
and gyroscopes on ground had not been applied. We 
present results obtained from XIS-0, one of the three front- 
side illuminated CCDs, because the XIS-0 has the best 
PSF suited for our purpose. The drifts of all four XIS im- 
ages due to the drifts of the satellite attitude were found 
to be almost identical. The results from XIS-0 are equally 
applicable to the other XIS data. 

3. Analysis 

3.1. Position Drift of Suzaku XIS Image 

In order to track a drift of a source position in 
the sky coordinates (X/Y), we determined the quasi- 
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AX det = AX cos 9 + AY sin 6 
AY dot = -AX sin0 + Ay cos 0, 



instantaneous peak of the sky image for every 200 second 
time bin by a two-dimensional Lorentzian fitting. The 
deviation of the peak position from the expected position 
(AX, AY) was computed at each time bin. The expected 
pixel position is (768.5, 768.5) if the source was aimed to 
be placed at the XIS nominal position. They were then 
converted into the detector coordinate system, (AXd e t, 
AYdet), using the following transformation: 

(1) 
(2) 

where 8 is the average roll angle, PA_NOM, stored in 
the event file headers. A similar method was adopted 
in Gotthelf ct al. (2000) for correcting the ASCA position 
accuracy. 

To illustrate the characteristics of the drift of the XIS 
images, the time history of a source position obtained from 
data of Her X-l (seq# 100035010) is shown in figure 2. It 
can be seen that the source position oscillates with am- 
plitudes of 10"-30" in both AXdet and AYd e t synchro- 
nized with the satellite orbital motion of a 96-min period. 
Furthermore, the center of oscillation has an offset toward 
the — AXdot direction by ~ 40", which causes a large as- 
trometric error in the Suzaku XRT-XIS system. 

We found that the position drift of the Suzaku XIS 
image can be well described by a wobbling motion syn- 
chronized with the satellite orbital motion and an offset 
which changes slowly in a much longer timescale than the 
orbital period during a fixed pointing observation. The 
offset appears to differ from one pointing to another. We 
derived the offset in both the X dc t and Ydct directions 
for each pointing by fitting the time-integrated XIS image 
with a two-dimensional Gaussian function. We denote 
the derived offsets as (AXdet) an d (AYdct), though these 
quantities are not necessarily the average of the above de- 
viations AXdet and AYdct- The offsets would be a good 
measure of the center of oscillation. 

In the left panel of figure 3 we show a scatter plot of 
the offset, (AXdot) versus (AYdct), derived for all of our 
samples. This plot represents the position accuracy of 
the Suzaku XIS image processed by the pipeline VI. 2. 
We found that the dispersion of the offsets is larger in 
the DETX (X det ) direction, ~ 70". The dispersion of the 
offset in the DETY direction is relatively small and the 
distribution tends to be positive. 

3.2. Parametrization of Drift 

We searched for related parameters to the offset and 
the wobbling motion of the X-ray image to recover the 
image that should be obtained in a condition free from the 
drift of the satellite attitude. We here assume that the 
long-timcscale offset is constant during a fixed pointing 
observation. 

3.2.1. Long-times cale offset 

First, we found that the offset in DETX, (AXdet), has 
a clear correlation with the ecliptic latitude of the source, 
f3 ec \. In figure 4 we present the distribution of (AXdet) of 
each observation against (3 ec \. A linear relation between 
the DETX offset (AXdet) and the ecliptic latitude [3 CC \ is 
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Fig. 2. (Top) Time sequence of AX^ ct (filled circles) and 
AYjjet (open circles) in the case of Her X-l. One pixel cor- 
responds to ~ 1.04". The origin of the source position is 
its nominal position, at which a source is expected to fall. 
(Bottom) Same as the top panel, but after the attitude-error 
correction (see subsection 3.2). 

clearly seen. The ecliptic latitude /3 ec i represents the an- 
gle of the XRT pointing direction from the ecliptic plane, 
and characterizes the way of illumination by the bright 
Earth. The strong correlation suggests that the geome- 
try between the Sun, Earth, and the satellite, plays an 
essential roll in determining the long-timcscale offset of 
the XRT pointing in the DETX direction. The config- 
uration of side panel #7 relative to the ecliptic plane is 
considered to somehow control the thermal distortion of 
the cross frames. 

The following linear relation is found to reasonably fit 
these data: 



(AX de t)(/3 ec i) = -38.5 



ecl 



-9.25 pixels. (3) 



^60 deg. 

This relation gives a DETX offset for a given observation 
using the ecliptic latitude of the target. It can be seen 
that (AXdet) — for f3 ec i ~ —14°. This is simply because 
the nominal position has been calibrated by MCG— 6— 30- 
15 (seq# 100004010), which has an ecliptic latitude of 
(3 ccl ~ -22?5. 
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Fig. 3. (Left) Offsets of XIS source positions from the nominal positions: (A-Xjet) versus (AY^et)- 1 pixel ~ 1.04". Aim points of 
the observations are placed at the XIS nominal position (filled circles) or at the HXD nominal position (open circles). (Right) The 
same plot after the position correction. 
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Fig. 4. Correlation between (iX^) and the ecliptic lati- 
tude /3 cc i of the sources. 

3.2.2. Wobbling with satellite orbital period 

We next explored the wobbling motion synchronized 
with the orbital motion of the satellite. After inspecting 
correlations with various parameters stored in the satellite 
housekeeping data, we arrived at two parameters related 
with the wobbling motion. One is a difference of temper- 
ature between radiators at the side panels #8 and #6, 
T_86 (in units of Kelvin), defined as 

T_86 = HK_XIS_RAD8_T1_CAL - HK_XIS_RAD6_T1_CAL, (4) 

where HKJCIS_RAD8_T1_CAL and HK_XIS_RAD6_T1_CAL are 

temperatures of radiators at the side panel #8 and #6, 
respectively. The newly introduced parameter is found to 
be a good indicator of the wobbling motion in the DETX 
direction. The temperature difference T_86 represents a 



temperature gradient in side panel #7 in the DETX di- 
rection. The other parameter is the elapsed time after a 
night-day transition, T_DY (in units of second), defined by 

f TJDY.NT for day-time, 

\ — TN_DY_NT for night-time, [ ' 

where day-time means a period in which Sun light irradi- 
ates the satellite. We make use of the parameter T_DY to 
parameterize the wobbling motion in the DETY direction 
because it is supposed to be coupled with thermal cycling 
of the orbiting spacecraft. 

Figure 5 shows a correlation between T_86 and AX^et, 
after removing the offset of (AXdct)- In the plot, we 
utilized only bright sources with a flux higher than 
2 counts s _1 among the samples in table 1, because the 
source position has to be determined for each 200-sec bin. 
For a fixed value of T_86 in a given observation, AXdct 
was found to have a scatter of about 10 pixels. We aver- 
aged AXdct for each temperature bin over each observa- 
tion. There is a tendency that the source position devi- 
ates from the nominal position to the +DETX direction 
as the temperature difference T_86 increases. We therefore 
parametrized AX^ct (pix) by the combination of a linear 
function of T_86 and (AXdet)(ftci) in equation (3) as 

AX dct = (AX dc t)(/3cci) + 0.4 x T.86. (6) 

In figure 6 we demonstrate a correlation seen between 
T_DY and AYdct- It appears that AYdet increases as T_DY 
increase for T_DY < 2000 and then starts to fall. We ap- 
proximated the relation between AYdct (pi x ) and T_DY 
with a broken line: 

( 10 + 0.007 X T_DY for T_DY < 0, 
AYdct = < 10 + 0.0025 x T_DY for < T_DY < 2000, (7) 
[ 45 0.015 x T_DY for 2000 < T_DY. 
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These empirical relations give an instantaneous devia- 
tion from the nominal position at each event time for a 
given observation, thus allowing corrections in the DETY 
direction. Note however that the absence of a direct mea- 
surement of the cross frame temperatures makes the cor- 
rections difficult. 
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Fig. 5. Correlation between T_86 and AX^ct — (AXj e t), de- 
rived by tracking the wobbling motion of bright sources. The 
filled points correspond to the XIS nominal cases and the 
open circles to the HXD nominal. The solid line represents 
equation (6), and the dashed lines show a range of scatter in 
(AXdct). 
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Fig. 6. Correlation between T_DY and AY dct . The filled 
points correspond to the XIS nominal cases and the open 
circles to the HXD nominal. The function defined by equa- 
tion (7) is drawn. 



4. Software Implementation 

We have developed a program, aeattcor, to correct the 
parameters of the satellite attitude using the parameteri- 
zation of drift of the XIS image described above. 

The program first calculates the expected deviation of 
the XRT pointing direction from that aimed at each event 



time using equations (6) and (7). The deviation corre- 
sponds to the error of the satellite attitude caused by 
any misalignment between the XRT-XIS system and the 
attitude control system. The program calculates Eulcr 
angles of the satellite attitude for the XRT-XIS system 
from those monitored on the attitude control system, and 
then updates the attitude file in which the Euler angles 
are stored. If one simply converts detector coordinates 
(DETX, DETY) to sky coordinates (X,Y) using the up- 
dated attitude file, the corrected sky coordinates (X,Y) 
can be obtained. A detailed description of the conversion 
scheme is given by Ishisaki et al. (2007). 

The program is now included in FTOOLS, the offi- 
cial analysis software released form the Suzaku Guest 
Observer Facility, and will be applied to the XIS event 
files released through the pipeline processing V2.0 or the 
later versions. 

5. Demonstration 

In this section we demonstrate the developed program, 
aeattcor, in application to the XIS images of sample 
sources. We show how degree the program improves the 
PSF and the position accuracy. 

5.1. Restored Image and PSF 

The left panel of figure 7 shows time-integrated XIS-0 
image of Her X-l before the attitude-error correction, as 
extracted from the VI. 2 event files. The image before the 
attitude-error correction looks as if there are two point- 
like sources separated by 30". It is because the position of 
Her X-l is oscillated with the satellite orbital period, and 
the source visibility is limited in selected orbital phases 
for Earth occultation (see figure 2). 

The right panel of figure 7 is an image after the attitude- 
error correction. One can realize that the deviation of 
the peak from the nominal position is improved, and the 
double-peak profile of the image before the correction dis- 
appears in the corrected image. In figure 2, we present a 
time sequence of the source position after the correction, 
where the large amplitude in the DETY is reduced. 

To quantify the improvement of the PSF core, we intro- 
duce a new parameter, "core sharpness" S, defined as the 
ratio of the averaged count per pixel within an inner circle 
of 10-pix radius to that within an outer annulus between 
10-pix and 20-pix radii. Here we determined the source 
centroid by fitting with a 2-dimcntional Gaussian func- 
tion. If the PSF is perfectly flat within 20 pixels from the 
center, the core sharpness is 1.0. For the instantaneous 
PSF free from the time-variable pointing direction, we ex- 
pect S ~ 2. We note that the "core sharpness" depends on 
each XIS because each XRT has a slightly different PSF. 
We here present the results of XIS-0. 

Figure 8 shows the distributions of S for the sources 
listed in table 1, before and after the attitude-error 
correction. The sources with a flux higher than 
50 counts exposure -1 chip -1 are shown separately with 
red histograms, because those bright sources are expected 
to severely suffer from pileup effects. It is clear that the 
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Fig. 7. Time-integrated XIS-0 images made with (X, Y) 
columns of Her X-l are displayed before (left) and after (right) 
the position correction. Counts per pixel normalized at the 
peak are shown on a linear color scale. North (+Y) is up and 
west (+X) is to the right. The double circles with a radius 
of 10" and 20" show the nominal (expected) position of the 
source. The source peak in the corrected image in the right 
panel is closer to the nominal position. Furthermore, a fake 
double-core of the original image is clearly converged. 

PSFs of the corrected image are sharpened and the core 
sharpness becomes close to 2, the value expected for the 
PSF free from the attitude error. 
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Core Sharpness, S 

Fig. 8. Distributions of core sharpness S for XIS-0. The top 
panel shows the distribution of S before the correction, while 
the bottom panel shows that after it. 5 = 1 corresponds to a 
flat PSF in the central 20-pix and S = 2 to the undistorted 
PSF core. Results for bright sources exceeding 50 counts 
exposure -1 , which are likely to be affected by photon pile-up, 
are separately shown in red. 

5.2. Improved Position Accuracy 

We determined the offset of the source centroid from 
the expected position for all sample sources before/after 
the attitude-error correction, in order to evaluate the im- 



provement on the source position accuracy. The distribu- 
tions of these offsets ((AXdet), (Aldet)) before/after the 
attitude-error correction arc plotted separately in figure 3. 
The improvement of the position accuracy is clearly seen 
from the difference of the two distributions. 

To quantify the position accuracy, figure 9 shows the 
distributions of Ar = yj (AATdct) 2 + (Aldct) 2 , i-e., the dis- 
tance between the source centroid and the expected po- 
sition, separately before/after the attitude-error correc- 
tion. Without the attitude-error correction, the devia- 
tion from the nominal position has a 90% error circle of 
Ar ~ 50" radius. The aeattcor restores the pointing ac- 
curacy to Ar < 19" (a 90% error circle). The improvement 
is brought mainly by the successful parameterization of 
the large DETX offsets. The residual position error after 
the attitude-error correction is somewhat worse compared 
with the preflight specification, ~ 10", and the position 
accuracy of ASCA, ~ 12" (Gotthelf et al. 2000). Given 
a lack of direct measurements of the temperatures of the 
cross frames, which are the most relevant parameters, we 
do not expect that further improvements can be easily 
made. We will monitor the Suzaku pointing accuracy to 
check its validity over the course of the mission. 
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Distance from Nominal Position (pix) 

Fig. 9. Distributions of the source distance Ar from the 
nominal position. The top panel shows the Ar distribution 
before the attitude-error correction, while the bottom panel 
shows that after correction. 

6. Summary 

We have developed an empirical method to restore the 
position accuracy and the PSF in the observed Suzaku 
XIS image. We analyzed 27 pointing observations of 
point-like sources using the XIS-0 data of the pipeline 
processing VI. 2 to investigate the characteristics of 
the drift of the XIS image which degrades the position 
accuracy and the PSF. We found that the deviation of the 
source centroid from the expected position, averaged over 
one pointing observation, can be well characterized by the 
relation between the deviation in the DETX direction and 
the ecliptic latitude of the pointing target. It is also found 
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that the oscillation of the source position synchronized 
with the satellite orbital period can be parameterized to 
some extent by the temperatures of the onboard radiators 
and the elapsed time since the night-day transition of 
the satellite. Based on the parameterization of the drifts 
of the XIS image, we have implemented new software, 
aeattcor, to correct the XRT pointing direction with 
respect to the attitude control system at each event time. 
By applying the attitude-error corrections to the XIS 
data of our sample, we demonstrated that the distorted 
core of PSF can be sharpened and showed that a 90% 
error circle of the sky position of a point-like source 
is reduced to 19" in radius. The results presented in 
this paper describe the current accuracy of performing 
astrometry with the Suzaku XIS. 

We thank the anonymous referee for useful comments 
and suggestions, which improved the paper. We are 
very grateful to Takashi Usui, Akira Okamoto, Kazuyo 
Mizushima, Tetsu Saitoh, Kazunori Shouji, Takayuki 
Tohma (NEC Co.) who gave us valuable comments on 
the thermal distortions of the Suzaku satellite. 
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Table 1. List of Suzaku Point Sources Used in This Study 
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Notes. — Units of right ascension are 


hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. AI cor 



and AY cor are pixel deviations of the source centroid from the nominal position in the corrected X/Y image, while AX and AY are for 
the uncorrected image. 

* Aim point of the observation. XIS: the target is placed at the nominal optical axis of the XIS. HXD: the target is placed at the nominal 
optical axis of the HXD, which is shifted in DETX by — 3.'5 from the XIS optical axis, 
t Roll angle in the observation, stored as PAJJ0M in the header of the event FITS files. 

t These sources have > 50 counts exposure -1 (0.4-10 keV) and consequently show photon pile- up in the central regions (within about 20 
pixels from the center) of the PSF. 



